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A determination of the excited energy eigenstates of the nucleon, s = I = ^, N , is presented 
in full QCD using 2 + 1 flavor PACS-CS gauge configurations. The correlation-matrix method is 
used and is built using standard nucleon interpolators employing smearings at the fermion sources 
and sinks. We develop and demonstrate a new technique that allows the eigenvectors obtained to 
be utilized to track the propagation of the intrinsic nature of energy-states from one quark mass to 
the next. This approach is particularly useful for larger dimension correlation matrices where more 
near-degenerate energy-states can appear in the spectrum. 



PACS numbers: 11.15.Ha,12.38.Gc,12.38.-t 



I. INTRODUCTION 

Resonances represent some of the rich dynamics of one 
of the fundamental interactions of Nature, the strong in- 
teraction of quarks and gluons. Lattice QCD is the only 
ab initio first principles approach to the fundamental 
quantum field theory governing the properties of hadrons 
and ultimately we wish to test our theoretical under- 
standing of resonances against their experimentally de- 
termined properties. 

From lattice QCD, the ground-state hadron spectrum 
is now relatively well understood However, gain- 
ing knowledge of the excited-state spectrum on the lat- 
tice presents additional challenges, as the excited energy- 
states are extracted from the sub-leading exponentials of 
the correlation functions. A determination of the excited 
state energy spectrum requires significant effort and some 
progress is now being made. We can expect the interplay 
between lattice QCD predictions and experimental mea- 
surement to be very productive in the coming years. 

In the case of nucleon resonances, the first positive 
parity excitation of the nucleon, the Ni~'^(1440) Pn or 
Roper resonance, has been a subject of considerable in- 
terest since its discovery in 1964 through a partial-wave 
analysis of pion- nucleon scattering data [Ij. This state 
has a surprisingly low mass, which is well below the first 
negative parity excitation. In constituent quark models 
with a harmonic oscillator potential this Pn state (with 
principal quantum number N = 2) appears above the 
lowest-lying odd-parity Sn (1535) state 0,13], whereas in 
Nature the Roper resonance is almost 100 MeV below the 
Sii state. This presents a phenomenological challenge 
to our understanding of level ordering. Similar difficul- 
ties occur with the = A* (1600) and i'^S;*(1690) 
resonances. Due to its surprisingly low mass, the Pn 
state has lead to enormous curiosity and much specula- 
tion about its nature. For example, the Roper resonance 
has been described as a hybrid baryon with explicitly ex- 
cited gluon field configurations [1, IB| , or as a breathing 
mode of the ground state Q or as a five quark (meson- 



baryon) state [8| . Significant resources have been devoted 
in the past from the lattice QCD perspective to find the 
elusive low-lying Roper state, in both quenched [9l-[25j 
and in full [26^301 QCD. 

The 'Variational method' [Sllls^ is the state-of-the-art 
approach for determining the excited state hadron spec- 
trum. It is based on the creation of a matrix of correlation 
functions in which different superpositions of excited- 
state contributions are linearly combined to isolate the 
energy eigenstates. A low-lying Roper state was identi- 
fied with this method using a variety of source and sink 
smearings in constructing correlation matrices [22l . [23 | 
in quenched QCD. Recent developments of algorithms 
and computational power have enabled the extension to 
full QCD. Some full QCD analyses using the variational 
method can be found in Refs. Here we 

consider the techniques of Refs. j22l |24| to explore the 
low-lying even- and odd- parity states of the nucleon us- 
ing 2+1-fiavor dynamical QCD gauge-field configurations 
from the PACS-CS collaboration [401 . A small subset of 
the results presented here have appeared in Refs. [29ll4ll|. 

The number of energy states revealed in the correla- 
tion matrix method depends on the number of unique 
operators chosen which have the quantum numbers of 
the desired states. Hence, a clear identification of these 
states is necessary to observe changes in these energy- 
states as a function of quark mass, in principle to the 
physical quark mass. This allows the quark mass depen- 
dence and structure of the extracted energy eigenstates 
to be explored systematically. The new technique that 
we develop here can be used when any parameter of the 
theory is varied to explore how the nature of the states 
and energies change with that parameter. 

The principal focus of this paper is to present the de- 
tails of our eigenvector analysis to track the states from 
the heavy to the light quark mass region. In doing so, 
we consider the states to illustrate the utility of the 
method. The operator basis is increased with the use 
of fermion source and sink smearings as in Ref. [29j . 
Then, the propagation of the energy states are pre- 
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sented after analyzing the state of eigenvectors at ad- 
jacent quarks masses. Results are presented for both the 
non-symmetric and symmetric eigenvalue equations, and 
these are compared and related to each other. 

If wc denote the dimensionality of the Hilbcrt space of 
the lattice Hamiltonian to be A^, then in an ideal world 
we would select N linearly independent operators to con- 
struct our N X N correlation matrix with sufficient sta- 
tistical accuracy and then diagonalize this to obtain the 
exact N energy eigenstates for this lattice Hamiltonian. 
Obviously and unfortunately this is not computationally 
feasible on any realistic lattice and the best that we can 
do is choose a relatively small number of operators, M, 
where M << N. If we choose these M operator interpo- 
lating fields wisely, then the subspace they span will have 
good overlap with the subspace spanned by the M lowest 
eigenstates of the Hamiltonian. If this is the case then 
with sufiicicnt statistics we can hope to extract good es- 
timates of the AI lowest energy states. If we observe, for 
example, that adding additional operators to increase M 
to M' reveals new low-lying states, then clearly we had 
not chosen our M operators wisely enough. The test of 
whether or not we have revealed all of the lowest states of 
the lattice Hamiltonian is that the process of adding ad- 
ditional and new operators docs not reveal new low-lying 
excited states. We present a clearer and more complete 
discussion of these issues in the Appendix. 

The paper is arranged as follows: Section |ll] contains a 
standard description of the mass extraction from a two- 
point correlation function with a brief introduction of 
Gaussian smearings at the fermion sources. The vari- 
ational method is presented in section IIIII followed by 
simulation details in section llVl Section IVl contains a dis- 
cussion of the energy eigenstates identification. Results 
for the flow of eigenvectors are presented in section IVII 
followed by concluding remarks in section IVIII Finally, 
a pedagogical discussion of the correlation matrix is pre- 
sented in Appendix in terms of the lattice Hamiltonian. 



II. ENERGY STATES FROM TWO-POINT 
CORRELATION FUNCTIONS 

A two point correlation function can be written as 

X 

where the Dirac spin indices are implicit. The operator 
Xj(0) creates states from the vacuum at space-time point 
and, following the evolution of the states in Euclidean 
time t, the states are destroyed by the operator Xii^) 
at the point {x,t). T indicates the time ordering of the 
operators. 

The energy eigenstates of hadrons are extracted using 
operators suitably chosen to have overlap with the de- 
sired states of interest. If we consider a baryon state B, 
then a complete set of momentum eigenstates provides, 



^ \B,p',s){B,p\s\^I, (2) 

B,p',s 

where B can also include multi-particle states that the 
operator x couples with. The substitution of Eq. ^ into 
Eq. (d]) yields 

{^l\x^{x)\B,p\s){B,p',s\xJm^)■ (3) 

Using the translational operator, the operator Xi{^) can 
be expressed as 

X^ix) = e^*e-^^--\.(0)e^^-^e-«*, (4) 

where H is the lattice Hamiltonian and P is the momen- 
tum operator whose eigenvalue is the total momentum p 
of the system. Inserting this into Eq. ([3]) we obtain 

X B.p'.s 

{n\x^mB,p',s){B,p',s\xJm^) 

= e-''^'Sp~j,,{n\x^mB,p',s){B,p',s\xJiOm) 

B,p',s 
B s 

The overlap of the interpolating fields x(0) and x(0) with 
positive and negative parity baryon states \B^) can be 
parametrized by a complex quantity called the coupling 
strength, Xg± , which can be defined for positive parity 
states by 



{n\xmB+,P,s) = Xg^. -^ubAp,s), (6) 

V ^B+ 



{B+,p,s\x{Om = Ab+W-^Sb+(p,s). (7) 
y J^B+ 

For the negative parity states one requires 

{n\x{0)\B- ,p,s) ^ Xs- ^l^l,us-{P,s), (8) 

{B-,p,s\xiOm) = -Ab-^/^Sb-(p,s)75. (9) 

Here, Xb± and A^i arc the couplings of the interpolat- 
ing functions at the sink and the source respectively and 
Mg± is the mass of the state B^ . Eb± is the energy 
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of the state B^, where Eg± = yA/^j_ Therefore, 

mass of a energy-state is obtained with the momentum 
projection of the correlation fmiction at p = 0. 

The standard spin sums may now be performed. For 
the positive parity hadron states, this can be expressed 
as 



1-P + Mb+ 



S 

and for the negative parity states, one encounters 



(10) 



^75 



^u§_(p,s)ug-(p,s) 



By substituting the above Eqs. for the positive and neg- 
ative parity states in Eq. ([5]) we obtain, 

Gij {t,p) = 2_^ Xb+ \b+ e B+ — 



E'B-XB-e--^-^ ^™-^^ ''-- . (12) 



the lower lying states. We employ a fixed boundary con- 
dition in the time direction for the fermions by setting 
Ut{x,Nt) = OVa? in the hopping terms of the fermion 
action with periodic boundary conditions imposed in the 
spatial directions. Gauge invariant Gaussian smearing 
[42] in the spatial dimensions is applied through an iter- 
ative process. The smearing procedure is: 



V'»(x,t) =^^^(f,x')V'»-i(^',i), 



(16) 



where, 



F{x,x') = (1 - a)(5x,x' + -^Y[Uf,ix)5^,^^+f, 

+ t/t(x-/i)<5,,. (17) 

where the parameter a ~ 0.7 is used in our calculation. 
After repeating the procedures A'sm times on a point 
source the resulting smeared fermion field is. 



At momentum p = 0, Eb± = Mb± , and a parity projec- 
tion operator T± can be introduced. 



r± = i(7o±i). 



(13) 



We can isolate the masses of the even and odd parity 
energy-states by taking the trace of Q with the opera- 
tors r_|_ and r_. The positive parity states propagate 
through the (1,1) and (2,2) elements of the Dirac ma- 
trix, whereas, negative parity states propagate through 
the (3,3) and (4,4) elements. The correlation function 
for positive and negative parity states can then be writ- 
ten as 



G±.(t,0)-Tr3p[r±e,,(t,0)] 



(14) 



The correlation function contains a superposition of 
energy-states, i.e. both ground and excited energy-states. 
The mass of the lowest energy-state, Mq± can be ex- 
tracted at large t where the contributions from all other 
excited-states are suppressed. 



\±\± p-M„±t 



(15) 



where, and are now couplings of interpolators to 
the lowest energy-state. 



A. Source Smearing 

The spatial fermion source-smearing (42j technique is 
applied to increase the overlap of the interpolators with 



III. VARIATIONAL METHOD 

The extraction of the ground state mass can be done 
straightforwardly using Eq. (fTS)) . However access to the 
excited state masses requires additional effort due to the 
presence of these energy-states at the sub-leading of the 
exponential. Here we consider the variational method 
[31I |32| . which allows for a variety of superpositions of 
excited-states in its cross-correlation discussed below. 

The variational method requires the cross correlation 
of operators so that the operator space can be diago- 
nalised and the excited state masses extracted from the 
exponential nature of the diagonalised basis. To access 
N states of the spectrum, one requires a minimum of N 
interpolators. With the assumption that only N states 
contribute significantly to Gij at time t, the parity pro- 
jected two point correlation function matrix for p = 
can be written as 

G±(t)^G±(i,0) 

= (X^Tr,p{r±(l]|x.(x)x,(0)|f^)}) 

X 

7V-1 

= E ^^^"6-"°*, (19) 



where Dirac indices are implicit. Here, Xf and A" are 
the couplings of interpolators Xi and Xj the sink and 
source respectively to eigenstates a = 0,-- - ,(A^— 1) and 
TOa is the mass of the energy-state a. The use of iden- 
tical source and sink interpolators provides A" = (A")* 
and then in the ensemble average Gfj{t) is a Hermitian 
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matrix, i.e. Gfj{t) = [G'^(i)]*. Moreover, considering 
both {u} and {u*} configurations makes Gfj(t) a real 
symmetric matrix. The N interpolators have the same 
quantum numbers and provide an A^-dimensional basis 
upon which to describe the states. Using this basis we 
aim to construct N independent interpolating source and 
sink fields which isolate N baryon states \Ba), i.e. 



N 



i=l 
N 



such that, 



(20) 
(21) 

(22) 
(23) 



where and z" are the coupling strengths of 0" and 0" 
to the state \Ba). Consider a real eigenvector which 
operates on the correlation matrix Gijit) from the right, 
one can obtain. 



G.AtK = (VTr,p{r±(f^|x,X,|f^)})«° 



For notational convenience, in the remainder of the dis- 
cussion the repeated indices i,j,k are to be understood 
as being summed over, whereas, a, which stands for a 
particular state, is not. 

In the ensemble average, Gij{t) = Gji{t). Therefore, 
considering ^[Gij{t) + Gji(t)] provides an improved un- 
biased estimator and enables the use a symmetric eigen- 
value equation as discussed below. To ensure that the 
matrix elements are all ^ 0(1), each element of Gij{t) 
is normalized by Gij{t)—^ (discussed in the 

\/Gii{0) y/Gj j{0) 

Appendix). 

In Eq. since the only t dependence comes from 

the exponential term, we can write a recurrence relation 
at time (to + At) as, 



G^jito + At)u" ^ e 



-m„ At 



Gy(toX, 



(25) 



for sufficiently large to and to + At [2ll. |43|. 

Multiplying the above equation by [Gij{to)]'^ from the 
left we get. 



[{G{to))-'G{to + At)]«" = 



— m^Ai a 

e " w 



(26) 



This is an eigenvalue equation for eigenvector u" with 



eigenvalue c" = e' 



-m„At 



We can also solve the left eigen- 



value equation to recover the eigenvector. 



<Gy(to + Ai) 



v?G,,{to) 



(27) 



Similarly, 



v"[G{to + At)iGito))-']=e- 



(28) 



Since Gij{t) is a real symmetric matrix v — u. The 
vectors and diagonalizc the correlation matrix at 
time to and to + At making the projected correlation 
matrix, 

vfG,j{t)uf^ = J^^z^z'^e"™"*. (29) 
The parity projected, eigenstatc projected correlator, 

<G±(i)<^GJ (30) 



is then used to obtain masses of different states. We 
construct the effective mass function 



M«ff(t) = In 



Glit,0) 
G%{t + 1,0)^ 



(31) 



and apply standard analysis techniques as described in 
Ref. [ill. 

Since G(to)"^^^G(to)^^^ 



: /, wc can rewrite Eq. ([26|) as 

Gito)-^Gito + At)G{toy^''^G{tof'^u'' = c"w" 

1 /2 

(24) Multiplying from the left by G{to) provides 



G(io)"'^'G(to + Ai)G(to)~'/'G(<o)'/'w" = c"G(to)'/'u" 
and defining, 

= G(to)^/\" (32) 

we find 

G{toy^'^-G{to + At)G{toy^'^w°' = c"w" (33) 
(also shown in Eq. (|A.8[) ). We note the matrix 

[G{to)-^'^G{to + At)G{to)-^'^] (34) 

is real symmetric, with the same eigenvalue c" and with 
the orthogonal to each other. If we had not used 
the [U + [/*] sum then the matrix in Eq. ([M]) would be 
hcrmitian and hence would still have real eigenvalues and 
orthogonal eigenvectors. The coefficients of interpolators 
creating an energy eigenstate is recovered by 



G{to) 



-1/2 



W 



(35) 



IV. SIMULATION PARAMETERS 



PACS-CS 2-1-1 flavor dynamical-fermion configura- 
tions [131 made available through the ILDG 44 1 arc used. 
These configurations use the non-perturbatively 0{a)- 
improved Wilson fermion action and the Iwasaki-gauge 
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action dH. The lattice volume is 32^ x 64, with (3 = 1.90 
providing a lattice spacing a = 0.0907 fm and lattice 
volume of « (2.90 fm)^. 

The degenerate up and down quark masses are 
considered, with the hopping parameter values of 
Kud = 0.13700, 0.13727, 0.13754, 0.13770 and 0.13781 
corresponding to pion masses of m,r = 0.702, 0.572, 
0.413, 0.293, 0.156 GeV [13; for the strange quark 
Ks = 0.13640. We consider an ensemble of 350 config- 
urations each for the four heavier quarks mass and 198 
configurations for the lightest quark. An ensemble of 750 
samples for the lightest quark mass is created by using 
well separated multiple fermion sources on each config- 
uration. We use the jackknife method to calculate the 
error, where the x^/dof for projected correlator fits is 
obtained via a covariancc matrix analysis. 

The complete set of local interpolating fields for the 
spin-^ nucleon are considered herein. Three different 
spin- flavor combinations of nucleon interpolators are con- 
sidered. 



X2{x) 



Xi{x) = e-''iu^''ix)Cj5d\x))u-{x) , (36) 

^abci-„.Ta 



e-={u''^^{x)Cd\x))-i^u'{x), 
{u^'' (x)C7574d'' {x))u'' {x) . 



(37) 
(38) 



The xi and X2 interpolators are used in Refs. jl2l.l46ll47|. 
The interpolator Xi is the time component of the local 
spin- 1 isospin-i interpolator which also couples to spin-^ 
states used, for instance, in Refs. fl6l [2lll48j. 

The local scalar-diquark nucleon interpolator, xi-i is 
well known to have a good overlap with the ground state 
of the nucleon. Also, this xi interpolator is able to ex- 
tract a low- lying Roper state in quenched QCD |2J| . On 
the other hand, the X2 interpolator has pseudoscalar- 
diquark structure in the nucleon, which vanishes in the 
non-relativistic limit, couples strongly to higher energy 
states. Each interpolator has a unique Dirac structure 
giving rise to different spin- flavor combinations. More- 
over, as each spinor has upper and lower components, 
with the lower components containing an implicit deriva- 
tive, different combinations of zero, one, two and three 
derivative interpolators are provided. 

The correlation matrices are constructed using differ- 
ent levels of gauge-invariant Gaussian smearing [l^ at 
the fermion sources and sinks [1^. A basis of smearing- 
sweep counts of 16, 35, 100 and 200 is selected following 
the extensive analysis of Ref. . 



V. EIGENSTATES IDENTIFICATION 

Let us consider M interpolating fields making an M x 
M parity-projected correlation matrix G{t). In solving 
the generalized eigenvalue equations of Eqs. ([^ and (^5]) 
we encounter the real and approximately symmetric ma- 
trices [{G{to))-^ G{to + At)] and [G{to + At) (G(to))"'], 
with the left and right eigenvectors and respec- 
tively. Thus the eigenvectors of these matrices are ex- 



pected to be approximately orthogonal (left table in Ta- 
ble As explained in the Appendix, the reason we have 
only approximate symmetry is that G{t) does not com- 
mute with itself at different times. This results because 
M < N. The more closely the subspace spanned by our 
M operators aligns with the subspace of the lowest M 
energy eigenstates of H, the less violation of symmetry 
there will be. If we do not use the [U + U*] sum, then all 
of the same arguments hold but with hermitian matrices. 
This feature enables the use of the generalised measure 

U'"^{mg,mg,) = u°'{mg)-u'^{mg,) (39) 

for the eigenvector , for example. This correlates eigen- 
vectors at different quark masses and may be useful in 
tracking states. 

In contrast, as already discussed, the matrix in Eq. ()34p 
is symmetric, hence the eigenvectors w°'(rng) are exactly 
orthogonal, i.e. w°'{mg) ■ ■w^{mg) = Sap (right table in 
Table m. 

Therefore, as in Eq. ([39|) . a generalised measure 



W^{mg,mg,)^w''img)-w^{mg,) (40) 

for the w" can be constructed to identify the states more 
reliably as we move from quark mass rrig to the adjacent 
quark mass m^' . 

In Table Hn the generalised measures U"^ {nig , m^' ) and 
I ) are presented for the heaviest and the sec- 
ond heaviest quark masses. It is evident that the off- 
diagonal elements of W"'^ (jUg , niqi ) are smaller than the 
U°'f^{mg,mg I ) and hence will be more reliable for identifi- 
cation and tracking of the energy eigenstates. Therefore 
we use yV"^(mg, nig') for this purpose. For each value of 
a there is only one value for /3 where the entry is within 
a few percent of 1. Thus this measure provides a clear 
identification of how eigenvectors in the hadron spectrum 
at quark mass are associated with eigenvectors at the 
next value of quark mass, mgi . 

Now we explain how we track the energy eigenstates 
from one quark mass to the next. Firstly, we label 
the extracted energy-states at the heaviest quark mass 
with a chosen set of symbols (most right column in 
Fig. [1]), where each symbol is assigned by the corre- 
sponding eigenvectors associated with it. These sym- 
bols are carried on to the lightest quark mass by look- 
ing at W"^ {mq,mgi) for adjacent quark masses going 
from the heaviest to the lightest. After the energy eigen- 
states are labeled at the heaviest quark, we look at the 
scalar product w"(mg) • w^{mgi) for the heaviest and 
the second heaviest quark mass (top left of Table IIIip . 
The scalar product shows that in this case all the diag- 
onal elements are larger than the off-diagonal, meaning 
there is no eigenvector crossing at these two heavier quark 
masses. A similar scalar product for the second heaviest 
to the third heaviest mass (top right of Table IIII[) shows 
that the fourth and fifth eigenvectors are crossed with the 
fifth and fourth at the third quark mass, and a similar 
situation for the sixth and the seventh. To illustrate our 
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TABLE I: The scalar product u"{mq) ■ u^{mq) (left) and 'w°'{mq) ■ w^{mq) (right), for the same quark mass, with 
four different levels of smearings. States are ordered from left to right and top to bottom in order of increasing 
excited-state mass, a and /? correspond to row and column, respectively. 



1.00 -0.18 0.02 -0.07 0.65 0.10 -0.32 -0.09 

-0.18 1.00 0.02 0.36 -0.10 -0.49 0.06 0.39 

0.02 0.02 1.00 0.15 0.07 0.06 0.42 0.03 

-0.07 0.36 0.15 1.00 -0.03 0.23 0.09 0.30 

0.65 -0.10 0.07 -0.03 1.00 0.15 -0.57 -0.13 

0.10 -0.49 0.06 0.23 0.15 1.00 -0.06 -0.61 

-0.32 0.06 0.42 0.09 -0.57 -0.06 1.00 0.17 

-0.09 0.39 0.03 0.30 -0.13 -0.61 0.17 1.00 



1.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 

0.00 1.00 0.00 0.00 0.00 0.00 0.00 0.00 

0.00 0.00 1.00 0.00 0.00 0.00 0.00 0.00 

0.00 0.00 0.00 1.00 0.00 0.00 0.00 0.00 

0.00 0.00 0.00 0.00 1.00 0.00 0.00 0.00 

0.00 0.00 0.00 0.00 0.00 1.00 0.00 0.00 

0.00 0.00 0.00 0.00 0.00 0.00 1.00 0.00 

0.00 0.00 0.00 0.00 0.00 0.00 0.00 1.00 



TABLE II: The scalar product u°'(mq) ■ u'^{mq') (left) and w°'{mq) ■ wl^ {mq,) (right), for k = 0.13700 and 

k! — 0.13727, with four different levels of smearings. States are ordered from left to right for m,/ and top to bottom 

for mq in order of increasing excited-state mass, a and j3 correspond to row and column, respectively. 



0.98 -0.29 -0.14 0.63 -0.07 0.10 -0.32 -0.08 

-0.19 -0.92 0.08 -0.03 0.14 0.06 0.42 0.05 

-0.16 0.07 0.99 -0.09 -0.04 -0.53 0.09 0.36 

0.63 -0.44 -0.02 0.99 -0.05 0.13 -0.55 -0.12 

-0.12 -0.11 0.40 0.00 0.75 0.00 0.08 0.36 

0.05 -0.11 -0.42 0.17 0.76 0.95 -0.12 -0.53 

-0.45 -0.17 0.03 -0.67 0.08 -0.05 1.00 0.18 

-0.09 0.00 0.34 -0.14 -0.34 -0.82 0.21 1.00 



1.00 -0.09 0.00 0.00 0.01 0.00 0.01 0.00 

0.09 0.99 -0.07 0.13 -0.01 0.00 0.01 0.00 

0.01 0.07 1.00 -0.01 0.00 -0.01 0.00 0.00 

-0.01 -0.13 0.02 0.98 -0.09 0.02 0.07 0.00 

0.01 0.01 0.00 -0.09 -0.97 0.21 -0.01 0.03 

0.00 0.00 0.01 0.00 0.20 0.95 -0.07 -0.23 

-0.01 0.00 0.00 -0.07 0.01 0.07 0.99 -0.01 

0.00 0.00 0.00 -0.01 -0.08 -0.21 0.01 -0.97 



analysis in our figures we track the eigenvectors from one 
quark mass to the next by connecting these similar eigen- 
vectors by lines. In Fig. [1] two lines (eigenvectors) cross 
at the second and the third quark mass. Wc then follow 
the above procedures for the third, fourth and fourth and 
fifth (the lightest) quark masses. 



It is well known in quantum mechanics the energies 
avoid level crossings as illustrated in Fig. [2l However, 
when two energy levels experience an avoided level cross- 
ing, the nature of the two eigenvectors is interchanged, as 
shown by the dotted lines in Fig. [21 In Fig.[Tl the first and 
second excited energy cigcnstates do not experience an 
avoided level crossing at pion mass of 413 MeV, whereas 
avoided level crossings are present for the third-fourth 
and the fifth-sixth excited energy-states. However, note 
that the avoided level crossings lie within the error bars. 
Results are presented as a function of quark mass (mq), 
with niqi = ArUq + mq where Am^ is small. In principle, 
as noted earlier a similar analysis can be performed for 
other lattice parameters in addition to the quark mass, 
such as the lattice spacing (a), volume {V), the lattice 
action etc. 



The eigenvectors are also tracked for the correlation 
matrix analysis with the xi ^^nd Xi interpolators and 
presented in Fig. [31 which can be compared with Fig. [H 



VI. QUARK-MASS FLOW OF EIGENSTATES 

A. Positive Parity 

A key feature of large correlation matrices is the ability 
to identify and isolate energy eigenstates which are nearly 
degenerate in energy. However, this approximate degen- 
eracy makes it difficult to trace the flow of states from 
one quark mass to the next. Thus a clear identification of 
these near-degenerate states through the features of the 
eigenvectors isolating the states is necessary in order 
to trace the propagation of the states from the heavy to 
the light quark-mass region. At this point it is useful to 
clarify our use of language. Where wc say eigenvector we 
arc referring to the orthogonal eigenvectors, w", of our 
symmetric (or hermitian) MxM correlation matrix G{t). 
Where we speak of energy eigenvalues and energy eigen- 
states, we are referring to the eigenvalues and eigenstates 
of the lattice Hamiltonian, H. Our goal in calculations 
is always to choose the M interpolators well so that the 
few (< M) lowest eigenvalues extracted are a good ap- 
proximation to the few lowest energy eigenvalues of H 
and so that M eigenvectors of our M x M correlation 
matrix G(t) capture the dominant characteristics of the 
corresponding eigenstate of H. 

The anticipated and relatively smooth fiow of the 
eigenvectors as a function of the quark mass is presented 
in Fig. m It appears that each eigenvector corresponds to 



7 



an energy eigenstate of H with the eigenvector cap- 
turing some of the core properties of the corresponding 
fuh energy eigenstate of H. While the quark-mass de- 
pendent trends can be significant, our approach reUably 
allows the identification of energy eigenstates at adjacent 
quark masses. 
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0.0 0.1 0.2 0.3 0.4 

m ^ (GeV^) 

FIG. 1; (Color online) N^'^ energy-states from 8x8 
correlation matrix of Xi,X2 interpolators from 
K = 0.13700 (to^ = 702 MeV, right-most column) to 
K = 0.13781 (to^ = 156 MeV, left-most column). The 
symbols follow the eigenvector as determined by 
considering the scalar product w°' ■ , as presented in 
Table IIIII Note that the dotted lines in the figure 
connect similar eigenvectors. Where these lines cross, of 
course the energy levels would not cross, but we would 
see an avoided level crossing as in Fig. [2] if we had data 
for every quark or pion mass. 



As the xi and X4 spin-flavor interpolators are very sim- 
ilar for the N^~^ channel, the overall flow of the eigen- 
vectors obtained from the Xi,X2 and X2,X4 correla- 
tion matrices are very similar in Fig. |4l Also, the over- 
all strength of the eigenvector components creating and 
annihilating N^~^ energy-states in the QCD vacuum re- 
mains approximately the same for the Xi;X2 and X2!X4 
cases (Fig. [5]), which implies that the eigenstate-energies 
isolated by the X1jX2 and X2jX4 analysis are the same. 
As the first excited state is purely xi-spin-fiavour domi- 
nated, this state is revealed in all the three different 8x8 
correlation matrix analyses. 

There are a few general trends apparent in Figs. El^a) 
and mb) which are worthy of note. Focusing on [D^a) 
for specific reference, we see that there is often compe- 
tition between different smearing levels in creating the 
states. A good example is in state two, where the 200 
sweep xi interpolator is complemented by the 35 sweep 
Xi interpolator at heavy quark masses, but the 35 sweep 
interpolator strength is phased out as one approaches 
light quark masses with strength transitioning to the 100 
sweep xi interpolator. Even in the ground state one can 




FIG. 2: (Color online) Illustration of an avoided level 
crossing. The solid lines illustrate how the energy levels 
avoid crossing, while the two dotted lines illustrate how 
the nature of the associated eigenvectors fiow. In the 
region where the energy states are closest each is an 
equal but orthogonal admixture of the two eigenvectors. 




0.0 0.1 0.2 0.3 0.4 

m ^ (GeV^) 



FIG. 3: (Color online) As in Fig. [H but with xi, X4 
interpolators. The dotted lines in the figure connect 
similar eigenvectors. 



see the importance of the 200 sweep interpolator increas- 
ing as one approaches the lighter masses. This effect is 
even stronger in Fig. [^b) for the X4X2 analysis where the 
100 sweep X4 operator is phased out in favour of the 200 
sweep X4- 

One can also observe the superposition of Gaussian 
smearings of different sizes being superimposed with rel- 
ative minus signs in a manner that will create nodes in 
the radial wave function of the interpolator of Eq. (|20p . 
Focusing again on Fig.[5]Ja) to provide a specific example, 
consider state 2. Here the widest Gaussian of 200 sweeps 
is complemented by a smaller Gaussian with the oppo- 
site sign. Moreover, as the quark masses become light, 
the smaller Gaussian grows in size. The result is that 
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TABLE III: The scalar product w"(m,) • w^{mq,), for k = 0.13700 (m^ = 702 MeV) and k' = 0.13727 (m^ = 572 
MeV) (top left), k = 0.13727 (m^ = 572 MeV) and k' = 0.13754 (m^ = 402 MeV) (top right), k = 0.13754 
(m^ = 402 MeV) and k' = 0.13770 (m^ = 293 MeV) (bottom left), k = 0.13770 (m^ = 293 MeV) and k' = 0.13781 
(to^ = 156 MeV) (bottom right), for an 8 x 8 correlation matrix of xi a-nd X2, with four different levels of 
smearings. States are ordered from left to right for m^/ and top to bottom for in order of increasing excited-state 
mass, a and /3 correspond to row and column, respectively. 



1.00 -0.09 0.00 0.00 0.01 0.00 0.01 0.00 

0.09 0.99 -0.07 0.13 -0.01 0.00 0.01 0.00 

0.01 0.07 1.00 -0.01 0.00 -0.01 0.00 0.00 

-0.01 -0.13 0.02 0.98 -0.09 0.02 0.07 0.00 

0.01 0.01 0.00 -0.09 -0.97 0.21 -0.01 0.03 

0.00 0.00 0.01 0.00 0.20 0.95 -0.07 -0.23 

-0.01 0.00 0.00 -0.07 0.01 0.07 0.99 -0.01 

0.00 0.00 0.00 -0.01 -0.08 -0.21 0.01 -0.97 

1.00 -0.04 -0.02 0.04 0.01 0.00 0.00 000" 

0.03 0.98 -0.21 0.04 -0.01 0.00 -0.03 0.00 

0.02 0.21 0.97 0.01 0.14 0.04 -0.02 -0.04 

0.01 -0.01 -0.13 -0.37 0.92 0.08 -0.03 -0.03 

-0.04 -0.04 -0.05 0.93 0.36 0.02 0.00 -0.01 

0.00 -0.03 -0.01 0.01 -0.01 -0.25 -0.97 -0.03 

0.00 -0.01 -0.04 0.01 -0.10 0.95 -0.24 -0.16 

0.00 -0.01 -0.02 0.00 -0.02 -0.15 0.07 -0.99 



1.00 -0.08 0.01 -0.01 0.01 0.01 0.00 0.00 

0.08 0.98 0.12 -0.03 0.09 0.01 0.00 0.00 

-0.02 -0.12 0.99 -0.08 0.00 0.00 0.00 -0.01 

-0.01 -0.09 -0.01 0.03 0.99 -0.10 0.00 0.00 

0.01 0.02 0.08 0.99 -0.02 0.01 0.07 0.05 

0.00 0.00 -0.01 -0.08 0.00 -0.08 0.99 0.07 

0.01 0.02 0.00 0.01 -0.10 -0.99 -0.08 0.03 

0.00 0.00 0.00 -0.04 0.00 0.03 -0.08 1.00 

1.00 -0.04 0.03 -0.02 0.01 0.00 0.01 0.00 

0.03 0.97 0.25 0.06 -0.02 -0.01 -0.01 -0.01 

-0.03 -0.24 0.94 -0.07 -0.21 0.00 -0.03 -0.02 

0.02 -0.06 -0.03 0.93 -0.36 -0.06 0.02 0.00 

-0.01 -0.05 0.20 0.35 0.89 -0.03 -0.21 -0.01 

-0.01 -0.01 0.06 0.04 0.18 -0.23 0.93 -0.20 

0.01 0.00 -0.02 -0.08 -0.04 -0.97 -0.22 0.05 

0.01 0.00 -0.04 -0.02 -0.04 0.00 -0.20 -0.98 



the radial node position of the wave function increases in 
distance as the quarks become lighter. 

Similarly, state 4 involves the superposition of three 
Gaussian smearings with alternating signs. Here the 200 
sweep interpolator is complemented with the 100 sweep 
interpolator with opposite sign which is complemented 
further by the 35 sweep operator, again with the opposite 
sign. Such a linear combination can create two nodes in 
the radial wave function. 

Finally, state 7 combines the 200, 100, 35 and 16 sweep 
interpolators with alternating signs such that a state with 
three nodes could be accessed. 

Turning our attention to the XiXi analysis, we see the 
flow of eigenvector components is not as smooth. Fig- 
ures Hl^c) and ^c) present the flow of the eigenvectors 
for this analysis. A careful comparison of the eigenstate 
spectrum with that from the X1X2 analysis at each quark 
mass and consideration of the eigenvector flow of the 
states reveals that the states dominated by xi in the 
X1X2 analysis are reproduced in the XiX4 analysis. How- 
ever, the remaining four states display a flow different 
from those revealed in the X1X2 analysis. Thus, a su- 
perposition of the xi,X2 and xiiX4 analysis provides 12 
unique energy states. 

In Fig. [HI a superposition of the two 8x8 analysis 
(8 X 8 X 2) of X1jX2 and Xi:X4 is presented. Scattering 
p-wave TriV and s-wave ttttN energy levels are also shown. 

For the two large quark masses, as seen in Fig. [6l 
the extracted lattice results sit close to the scattering 
two particle p-wave Ntt threshold {E^ + i5^) with back- 
to-back momenta, p = (27r/L2,, 0, 0) and s-wave Ntttt 
threshold (A/jy + A/^ + M^^), whereas the masses for the 
lighter three quarks sit much higher. There is no evidence 



of these scattering states at light quark masses. 

It is noted that the couplings to the multi-particle 
meson-baryon states are suppressed by 1/VV relative to 
states dominated by a single-particle state. Due to the 
large volume of our lattice, it is likely that multi-particle 
states will be suppressed and missed in our spectrum, 
particularly at lighter quark masses where the quark- 
mass effect also acts to suppress the spectral strength. 
Further analysis of finite volume effects [4§| on the spec- 
trum is highly desirable. Future calculations should also 
investigate the use of five-quark operators to ensure bet- 
ter overlap with the multi-particle states and to better 
resolve and probe the excited state spectrum [s^ [s^] . 



B. Negative Parity 

Now we can repeat our analysis for the identification of 
the iVi~ states. In Table HVl the scalar product W"'^ = 
w"(mg) • w^{mqr) for all the quark masses is presented. 

In Figs. [7] and El the spectrum from the 8x8 

analysis involving X1jX2 and X1jX4 is presented, respec- 
tively. While the XiiX4 analysis is able to extract a low- 
lying energy state, it misses the near-degenerate second 
energy state in this channel. This second energy-state 
is revealed in the xi j X2 spin-flavor combination present- 
ing two nearly-degenerate low-lying states, which is in 
accord with the quark model based on SU{6) spin-flavor 
symmetry. Recall that three spin-i quarks may provide 
a total spin of s = ^ or | , the L = 1 state can couple two 
different ways to provide a J ~ ^ state, hence providing 
two orthogonal spin-i states in the L = 1. 70 plet repre- 
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(a) Eigenvector components for an 8 X 8 correlation matrix 
with xii X2 interpolators. Odd and even numbers in the 
legend correspond to the xi and X2 respectively. 




-1.0 



EVect 1 EVect 2 EVect 3 EVect 4 EVect 5 EVect 6 EVect 7 EVect 8 

(a) Eigenvector components for an 8 X 8 correlation matrix 
with XI I X2 interpolators. Odd and even numbers in the 
legend correspond to the xi and X2 respectively. 
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(b) As in Fig. I4al but for the X2 and X4 interpolators, 
and even numbers in the legend correspond to the X4 
X2 respectively. 
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(b) As in Fig. I5al but for the X2 and X4 interpolators, 
and even numbers in the legend correspond to the X4 
X2 respectively. 
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(c) As in Fig. I4al but for the xi and X4 interpolators, 
and even numbers in the legend correspond to the xi 
X4 respectively. 
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FIG. 4: is presented for the five different quark 
masses for the N^'^ channel after identifying 
eigenvectors via w°'{'mq) ■ w^{mqr). For each eigenvector 
shown in the horizontal axis, the eigenvector 
components are plotted in order of increasing quark 
mass from left to right. Note that Evect 1 to Evect 8 
correspond to eigenvectors to w^. In the legend, 
subscripts (1, 2), (3, 4), (5, 6) and (7, 8) correspond to 
the smearing-sweep levels of 16, 35, 100 and 200, 
respectively. 
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(c) As in Fig. I5al but for the xi and X4 interpolators. Odd 
and even numbers in the legend correspond to the xi and 
X4 respectively. 

FIG. 5: u" is presented for the five different quark 
masses for the N^~^ channel. For each eigenvector 
shown in the horizontal axis, the eigenvector 
components are plotted in order of increasing quark 
mass from left to right. Note that Evect 1 to Evect 8 
correspond to eigenvectors to u^. In the legend, 
subscripts (1, 2), (3, 4), (5, 6) and (7, 8) correspond to 
the smearing-sweep levels of 16, 35, 100 and 200, 
respectively. 
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FIG. 6: (Color online) Nj^ energy-states from 
8x8x2 correlation matrices of Xi,X2;X4 [l^- The 
p-wave Ntt scattering threshold (with one unit of lattice 
momentum) and the s-wave Ntttt threshold are 
presented by dashed and dotted lines, respectively. 



sentation of SU{Q). Both of these states have a width of 
«150 MeV. 

As in Fig. |31 the eigenvector components for differ- 
ent quark masses are presented in Figs. [9l and [TOl for the 
channel. It is interesting to note that the scalar- 
diquark interpolator xi dominates the lowest two N~ 
states (Fig. [T0|) when available. The X2 interpolator 
makes an important contribution in creating the second 
energy state in this channel. 

In Fig. [TTl a superposition of the two 8x8 analysis 
(8 X 8 X 2) of X1jX2 and Xi:X4 is presented. Scattering 
p-wave TriV and s-wave TTTriV states are also shown. The 
results for the lowest energy-state at the two heavier pion 
masses sit close to the scattering s-wave N+tt (M^+A^w) 
threshold indicating that these results may be scattering 
states at these pion masses. However, they disappear 
from our spectrum at the light pion masses. A similar 
situation also prevails the second energy-state, where the 
state sits close to the p-wave En+Ett+AItt and ETr+E-^ + 
Mjv scattering threshold with back-to-back momenta of 
one lattice unit, p~ (27r/_L^, 0, 0). 



VII. CONCLUSIONS 

In this paper, a comprehensive analysis for the nucleon 
spectrum with / = 1/2, s = 1/2, is presented using the 
correlation matrix approach. In particular, a method for 
encrgy-cigcnstate identification and flow is presented and 
demonstrated for the positive and negative parity chan- 
nel. Details of the method developed for an identifica- 
tion and the propagation of the energy-states from heavy 
to light quark mass region are provided. In particular, 
the new technique is useful in identifying the flow of the 
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FIG. 7: (Color online) Masses of energy-states 
from an 8 X 8 correlation matrix of Xij X2 in Table ITVl 
The dotted lines in the figure illustrates the eigenvector 
fiow. 




0.2 0.3 
m ^ (GeV^) 

FIG. 8: (Color online) As in Fig. [3 but with xi,X4 
interpolators. The dotted lines in the figure illustrates 
the eigenvector flow. 



near-degenerate energy eigenstates from one quark mass 
to the next. The eigenvectors obtained from the eigen- 
value equations for several 8x8 correlation matrices are 
utilized in tracking the eigcn-energy states. 

In presenting the results, both non-symmetric left and 
right eigenvalue equations and a symmetric eigenvalue 
equation are considered. While the masses are the same 
in the two different approaches, the eigenvectors obtained 
from the symmetric matrix are orthogonal. Thus the 
generalized measure W"^ is used to track the flow of 
eigenvectors with quark mass. The scalar product of the 
eigenvectors shows its robustness in tracking the flow of 
the energy eigenstates even when the energies are nearly 
degenerate. 

The coefficients of the interpolators creating and anni- 
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TABLE IV: The scalar product w°'{mq) ■ w^(mq'), for k = 0.13700 {m^ = 702 MeV) and k' = 0.13727 (m^ = 572 
MeV) (top left), n = 0.13727 {m^ = 572 MeV) and k' = 0.13754 {m^ = 402 MeV) (top right), k = 0.13754 
(m^ = 402 MeV) and k' = 0.13770 (m^ = 293 MeV) (bottom left), k = 0.13770 (m^ = 293 MeV) and k' = 0.13781 
(m^ = 156 MeV) (bottom right), for an 8 x 8 correlation matrix of xi a-nd X2, with four different levels of 
smearings, for the states. States are ordered from left to right for m^/ and top to bottom for niq in order of 

increasing excited-state mass, a and /? correspond to row and column, respectively. 



-0.03 0.99 -0.08 -0.04 0.05 -0.01 0.00 0.01 

1.00 0.02 0.01 -0.08 -0.02 -0.01 0.00 0.01 

0.00 0.07 0.97 0.05 0.22 -0.03 -0.04 0.00 

0.07 0.03 -0.07 0.96 0.01 -0.27 0.00 -0.02 

-0.02 0.06 0.20 0.02 -0.95 -0.01 -0.23 -0.01 

0.03 0.02 0.01 0.26 0.00 0.94 -0.01 -0.21 

-0.01 0.02 0.08 0.01 -0.22 0.01 0.97 0.01 

0.00 0.00 0.00 0.08 0.00 0.20 -0.02 0.98 

0.91 0.40 0.02 0.02 0.01 -0.05 0.00 0.00 

0.40 -0.91 0.00 0.01 -0.02 0.01 -0.01 0.00 

-0.01 -0.01 0.96 -0.27 0.01 -0.01 0.00 0.02 

-0.03 0.00 0.27 0.96 0.01 0.01 0.02 0.00 

0.04 0.03 0.01 -0.01 -0.22 0.97 0.02 0.01 

0.01 -0.01 -0.01 -0.01 0.98 0.22 0.04 0.00 

0.00 0.00 -0.02 0.01 0.01 -0.01 -0.12 0.99 

0.01 -0.01 0.00 -0.02 -0.04 -0.03 0.99 0.12 



-0.22 0.97 0.04 -0.10 -0.01 0.04 0.01 -0.02 

0.97 0.22 -0.07 -0.07 0.01 0.02 -0.02 -0.01 

0.07 -0.03 0.99 -0.04 0.07 0.00 -0.02 -0.01 

0.05 0.11 0.04 0.99 -0.01 0.05 0.01 -0.04 

-0.01 0.02 -0.07 0.02 0.99 -0.09 0.04 0.03 

0.02 0.04 0.01 0.04 -0.09 -0.99 -0.02 -0.02 

0.02 -0.01 0.01 -0.01 -0.04 -0.01 0.99 -0.13 

0.01 0.02 0.02 0.03 -0.03 -0.02 0.12 0.99 

0.98 0.17 -0.06 -0.01 0.01 0.01 0.00 0.00 

-0.17 0.99 -0.01 -0.03 0.01 0.01 -0.01 0.00 

0.05 0.03 0.77 0.64 -0.01 0.00 0.00 0.02 

0.04 -0.01 0.64 -0.77 -0.04 0.00 0.01 0.00 

0.00 0.00 0.02 -0.01 0.66 -0.75 -0.04 -0.01 

-0.01 -0.01 0.03 -0.02 0.75 0.66 -0.04 -0.02 

0.00 0.01 0.00 0.01 0.05 -0.01 0.96 -0.26 

0.00 0.00 -0.01 -0.01 0.04 0.01 0.26 0.96 



hilating a state in the QCD vacuum arc also presented. 
The flow of the eigenvectors reveals a smooth pattern and 
presents important insights into baryon structure and its 
evolution with quark mass. 

Another interesting result of this paper is that, the cor- 
relation matrix method can be used to track the energy- 
states that are involved in an avoided level crossing. It is 
noted that the avoided level crossings lie within the er- 
ror bars, but the demonstration of the robustness of the 
approach remains. 



\En) be a complete orthonormal basis of energy cigcn- 
statcs, then the state \xi) can be written as 

N N 
\X^) - E \E^)(E^\X^) = E^'J-I^^-)' 

where, = {Ej\xi), which is analogous to Eq. ([71). In 
matrix form Eq. (jA.ip can be written as 
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Appendix: Pedagogical Discussion of the Correlation 
Matrix 

1. For M ^ N 

Let us consider an N dimensional Hilbert space with 
a hamiltonian, H, and let \xi), \X2), ■ ■ ■ , \xn) be N lin- 
early independent states. Similarly, let \Ei), |i?2), 



" Ixi) " 




' \Ei) ' 




= [C] 




. \xn) _ 




_ \En) _ 



Since \xi) is linearly independent, then C must be non- 
singular and so there must exist {C~^)ij such that 



" \Ei) ' 




' \Xi) ' 












. \xn) _ 



Similarly, {xi\ can be expressed as {xi\ = 

J:U(^^\E,){E,\^J:U^^AE^\■ 

Let us now define \xi) = xll^) for N linearly indepen- 
dent field operators (xj), where the dagger denotes the 
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(a) Eigenvector components for an 8 X 8 correlation matrix 
with xii X2 interpolators. Odd and even numbers in the 
legend correspond to the xi and X2 respectively. 
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(a) Eigenvector components for an 8 X 8 correlation matrix 
with XI I X2 interpolators. Odd and even numbers in the 
legend correspond to the xi and X2 respectively. 
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(b) As in Fig. I9al but for X2 and X4 interpolators. Odd and 
even numbers in the legend correspond to the X4 and X2 
respectively. 
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(c) As in Fig. I9al but for xi and X4 interpolators. Odd and 
even numbers in the legend correspond to the xi and X4 
respectively. 
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(b) As in Fig. llOal but for X2 and X4 interpolators. Odd 
and even numbers in the legend correspond to the X4 and 
X2 respectively. 
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(c) As in Fig. llOal but for xi and X4 interpolators. Odd 
and even numbers in the legend correspond to the xi and 
X4 respectively. 



FIG. 9: w"" for the five different quark masses are 
presented for th.e channel after identifying 

eigenvectors via w°'{'mq) ■ w^{mqi). For each eigenvector 
shown in horizontal axis, the eigenvector components 
are plotted in order of increasing quark mass from left 
to right. Note that Evect 1 to Evect 8 correspond to 
eigenvectos to w^. In the legend, subscripts (1, 2), 
(3, 4), (5, 6) and (7, 8) correspond to the 
smearing-sweep levels of 16, 35, 100 and 200 
respectively. 



FIG. 10: for the five different quark masses are 
presented for the channel. For each eigenvector 

shown in horizontal axis, the eigenvector components 
are plotted in order of increasing quark mass from left 
to right. Note that Evect 1 to Evect 8 correspond to 
eigenvectos to u^. In the legend, subscripts (1, 2), 
(3, 4), (5, 6) and (7, 8) correspond to the 
smearing-sweep levels of 16, 35, 100 and 200 
respectively. 
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FIG. 11: (Color online) Masses of iV^ energy-states 
from 8x8x2 correlation matrices of xi, X25 X4- 

adjoint. Then we may write 

N 

k,l = l 
N 



k=l 
N 



kj 



k=l 



or, 



/ e--^i* 



m)] = [c* 



-Eit 





[C 



*ti 



V 



(A.4) 



where, [E{t)] is obviously diagonal. The Eq. (|A.4p is 
analogous to Eq. (fT9|) . From this point on for notational 
convenience we will no longer use square brackets to de- 
note matrices. We see that, G{t) = G{ty for all t. From 
Eq. (|A.4|) . we can also write 



(A.5) 



where, C* = C*y^E(Jo) and C"*t ee y^E(Jo)C*l Simi- 
larly, we can write 



G(to + At) = C'*E{At)C'*K 



(A.6) 



Let us consider the polar decomposition of C* as 
C* = \/C'*C"*t t/ = ^/GjtojU, where [/ is unitary 



Similarly, C'*'< = W ^/G{t^. Then, Eq. ([Xe]) can be 
written as 



G(to + At) = VGMt^^^(At);7tv/G(M (A.7) 
or equivalently 



G(At) = v/G(to) G(to + At)VG(to) 

= UE{At)U\ {A.. 

where we are using the notation 



E{At) 



p-E^At 
e-^^^* 











(A.9) 



Denoting the normalized eigenvectors of G{At) as 
for i = I,-- - ,A^, then U consists of columns U = 



-;ii„7t2i i„7?Afi Hence we see that 



G(At)w^ = e 



(A.IO) 



Multiplying Eq. ([AIQ]) by G(to)"^/^ from left and defin- 
ing 



gives 



G{to)-^G{to + At)?? = e-^-^V, 
which is analogous to Eq. 



(A.ll) 
(A.12) 



2. For M < N 

Let xij X2: • • • , be M linearly independent inter- 
polating field operators with M < N. Then we may 
consider these M Xi's as a subset of a complete set of N 
interpolating operators. Then as before, we can define 
G^^^'^it) as an M X M correlation matrix, then G^' {t) 
can be written as the upper left M x M block of G(t) 
such that 



G(t) 



G^^(t) G"(t) 
G'"(t)t G(^-^^)(t) 



(A.13) 



where the off-diagonal rectangular matrix G'"(t) has el- 
ements G™(t) for i = 1, • • • , M and j = M 1, • • • , A^. 
Clearly G™(t) mixes the upper and lower diagonal blocks. 

In terms of the full N x N correlation matrix G(t) we 
define G(At) as before and it has the form 



G(At) 



G^^(At) G„(At) 
G,„(At)t G(^-^^)(At) 
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which is of course diagonalizcd by the full N x N unitary 
matrix U as shown in Eq. (jA.Sp . 

Let us now temporarily assume that the ofF-diagonal 
elements of G{At) are zero, then 



nient to redefine 



G{At) 



G^^(Ai) 

G(^-^^)(At) 



and then it follows that the time independent unitary 
matrix can be written as 



U = 



U^^ 

U^N~M) 



This will occur if and only if we have chosen our M in- 
terpolating fields such that they span exactly the same 
subspace as M of the exact energy eigenstates. 

But we will certainly never exactly achieve this and so 
mixing will occur through G'"'(t) 7^ 0. In this case U will 
not have this convenient block diagonal form. Let U' di- 
agonalize G^^(At), i.e. U' is an AI x M unitary matrix 
defined at this particular At. Then if G"^{t) 7^ we see 
that U' and in general U' will not be indepen- 

dent of At. Note that M linearly independent interpo- 
lating field operators will operate on the vacuum and give 
rise to M linearly independent states. The more closely 
these M linearly independent states come to spanning 
the same subspace as M exact energy eigenstates, then 
the smaller will be G™'{t) and the more block diagonal 
will be G{At) and U. We can only ensure that G^^(At) 
commutes with itself at all At if U is block diagonal, i.e., 
if G™(t) = so that U' = and is therefore time 
independent. Then the smaller will be the ofF-diagonal 
elements and the smaller will be the mixing and the ex- 
cited state contamination. 

In practice M << N, since on a lattice the dimension- 
ality of Hilbert space is in the many millions. We would 
like to choose the M interpolating fields such that they 
span an M-dimensional subspace that has the greatest 
overlap possible with the subspace spanned by the M 
lowest energy eigenstates \Ei), • • • , \Ei\i). With the low- 
est energy eigenstates our numerical errors will be mini- 
mized, since by working at large Euclidean times we min- 
imize the influence of the higher excited state contami- 
nation and so optimize our extraction of the lowest M 
energy eigenstates. 

Since the interpolating fields produced states \xi) that 
were not normalized in any way, it is numerically conve- 



where there is no sum over i. Therefore, we attempt to 
put the strengths of our interpolating fields at a compa- 
rable level by defining 



1 



1 



vcm^'' yGf(o)' 

to ensure that the matrix elements of G'^^{t) are all 
~ 0(1) to maximize numerical significance of all 
combinations. This is completely legitimate as it is sim- 
ply a change to the normalization chosen for our interpo- 
lating fields, which we are free to do. The matrix G'^\t) 
is hermitian except for the effects of finite ensemble and 
round-off errors. 
Let us define 

G"{t) = ^(G'^(t) + G'^^\t)), (A.14) 

as an improved unbiased estimator of the ensemble aver- 
age for G'^{t). Then G^{t) is exactly hermitian. There- 
fore, as before, we may diagonalize 



d"'{At) = [G^^(to) '^^G''{to + A<)G^(to)" 



— 1/2 

which is also hermitian, as G^^to) and G*^(to + At) 
arc obviously hermitian. The eigenvectors {w") obtained 
from diagonalizing the above matrix are therefore or- 
thonormal. It is noted that using the U + U* trick p^ . 
where these J7's are links here, the hermitian correlation 
matrix G^\t) is real symmetric and so the eigenvalues 
remain real and the eigenvectors orthogonal. Again, it 
is important to note that the more poorly we choose 
our M interpolating fields then the bigger will be the 

off-diagonal elements of G(t). Hence the more time- 
dependent will be the U' and the less reliable our ex- 
tracted energies and eigenvectors. 

Now we can consider G^^ (t) as a function of quark 
mass, TOq, to identify how a given state evolves with 
iriq. If = TTiq + Auiq with ArUq small, we expect 

(mq) .w^ {m'y) ^ 1 a i = j and w'' {mq) .w^ (m'^) ^ if 
i 7^ j- 



(A.15) 
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